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In the thermodynamic limit two distinct states of matter cannot be analytic continuations of each other. 

Classical phase transitions are characterized by non-analyticities of the free energy. 1 1] For quantum 

\Q , phase transitions |2] (QPTs) the ground state energy often assumes the role of the free energy. But in 

f^ ' a number of important cases this criterion fails to predict a QPT, such as the three-dimensional metal- 

f^ ' insulator transition of non-interacting electrons in a random potential (Anderson localization) |3]. It 

0^^ . is therefore essential that we find alternative criteria that can track fundamental changes in the inter- 

» , ■ nal correlations of the ground state wavefunction. Here we propose that QPTs are generally accompa- 

^ l ' nied by non-analyticities of the von Neumann (entanglement) entropy. |4] In particular, the entropy is 

.^ , non-analytic at the Anderson transition 1 5], where it exhibits unusual fractal scaling. We also exam- 

I ine two dissipative quantum systems of considerable interest to the study of decoherence and find that 

^O ' non-analyticities occur if and only if the system undergoes a QPT. [6] 

1— H ' n n n i— i 

^ • Recently there has been intense interest in using entanglement to study quantum phase transitions. |7, 8, 9, 10] Since entan- 

s" ' glement is a unique quantum phenomenon which characterizes wave functions that cannot be factored into single particle states, 
"^ one expects that it should play a role in generating the correlations that exist on all length scales at a quantum critical point. But 

. ; a state can be entangled without being critical — consider, for instance, the singlet state of two spin-1/2 particles. So it is not 

d mere entanglement which concerns us here, but also the criticality of infinitely many interacting degrees of freedom. 

H In an important paper on black hole entropy, Bekenstein 1 1 1 ] demonstrated the usefulness of information entropy. The concept 

is readily applicable to a quantum mechanical ground state. If we consider the the entire wave function, the conventional entropy 
is of course zero (assuming the state is non-degenerate), and all we can study is the analyticity of the ground state energy. In 

Q ■ those cases where this is a smooth analytic function of a tuning parameter that drives a QPT, we get no useful information. 

C_) ' Yet in many cases we know that the wave function represents special internal correlations, as in Laughlin's quantum Hall wave 
function, LlZl or in the Bijl-Dingle-Jastrow-Feynman wave function LI 3.1 describing the superfluid state of bosons. How can we 
,__! quantify such correlations, in particular when they change across a transition? 

J> Let us divide a system into two pieces and calculate the entanglement between them. In the ground state | $) , the density matrix 

04 of the full system is p = | ^I*) (vj/ 1 . Now partition the system into two parts A and B, where A denotes the subsystem of interest and 
'^ . B the environment whose details are of no interest. We construct the reduced density matrix pA by tracing over the environmental 
,,j. '■ degrees of freedom, akin to integrating out the microstates consistent with a given set of macroscopic thermodynamic variables. 
^\ • The von Neumann entropy, defined by S* = — Tr(/9A In Pa), provides a measure of bipartite entanglement and therefore contains 
\Q ' information about the quantum correlations present in the ground state. Note that S = —Tt{pa In p^i) = —Tt{pb In pb), where 
^^ > the reduced density matrix pB is obtained by tracing over the degrees of freedom in subsystem A. 

We consider three important models. In addition to demonstrating the use of von Neumann entropy, these models illuminate 
the role of disorder and dissipation in entanglement, inter alia in quantum computation. The two dissipative systems that we 

Ch discuss have been recently addressed from this point of view, [14] but the results are incorrect. Our first example, however, is 

F^ . Anderson localization in three dimensions (d ~ 3), which has been studied extensively and is known to have a quantum critical 

P^ ' point. A critical value of disorder separates a system where all states are localized from one where some states are extended. At 

O the critical point the wave function exhibits a fractal character 1 15] 

J-J The relevant Hamiltonian for a d-dimensional lattice, expressed in terms of fermionic creation (c]) and annihilation (ci) 

^ . operators, is — suppressing spin indices, irrelevant in this case — 

^ — — 

where the sum in the second term is over distinct nearest-neighbor pairs, and the site energies Ei are chosen from a uniform 
random distribution between —14^/2 and W/2. In what follows, we shall choose the unit of energy such that t = 1. Since this is 
a one-electron problem, the single-site density matrix in the state a is 

p? = icP(|i>(i|). + (i-iV'rP)(io)(oi). (2) 

where ip'j' is the amplitude of an energy eigenstate a at site i. The ket |0) (|1)) corresponds to the site being unoccupied 



a 



if = ^i?4c,-t^(ctc,+h.c.), (1) 



(occupied). The von Neumann entropy, averaged over all the eigenstates A and the number of sites N, is 

^ N A 



AN 



(3) 



n— 1 a— 1 



For a random system S must be further averaged over realizations of the random potentials. The resulting quantity, denoted by 
S, measures the correlation between the amplitude at a given site and the rest — local entanglement, llldl 
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FIG. 1 : The data collapse of the calculated von Neumann entropy of the d = 3 Anderson localization problem. The exponent v and Wc were 
set to the known values 1.35 and 16.5, respectively. The exponent y was then determined to be 2.6. The system sizes are indicated in the inset. 
The systems of sizes 7 and 9 were averaged over 10 different realizations of disorder, while the system of size 11 was averaged over 5 
realizations. For the 13^ lattice, it did not seem necessary to average over multiple realizations. 



We compute S by diagonalizing the Hamiltonian for finite-sized systems in d = 3 and then attempt to collapse the data using 
the finite size scaling hypothesis for the singular part of the entropy: 11711 



L-yML'/'-w), 



(4) 



where L — N^^"^ is the length of the system, w ^ \W — Wc\/Wc, and Wc is the critical disorder at which Anderson localization 
occurs. The critical exponent i/ is the exponent for the localization length, the value of which is known to be approximately 1.35; 
Wc is known to be approximately 16.5. 1 18|The universal scalin g fu nctions f± refer to W greater or less than Wc- It has been 
conjectured that for dimensions d > 2, the exponent y = (d — 1). lIlTII Examining the data collapse in Fig. [2 we find instead that 
the exponent y — 2.6, cl fractal dimension less than three. This is similar to the fractal scaling of the participation ratio. 1 15l ll9ll 
The origin of the fractal dimension is the fractal nature of the wave function at the mobility edge separating the delocalized states 
from the localized states. Of course, data collapse without an objective criterion can be treacherous. An interesting criterion 
used to justify the goodness of collapse is shown in Fig. |2]. The quantitative measure is obtained from the minimization of a 
specially defined function of the residuals called Pb, |20j which satisfies Pb > Pb|absmin = 0. 

The fractal dimension of the von Neumann entropy can be confirmed from another measure of entanglement called the 
linear entropy, Il21ll which is 5^ = 1 — Trp^. For the present problem, it is easy to show that the participation ratio P^^-' ~ 
\ Tlii a IV'f 1^ = 1 ~ Sl/2. This is to be averaged over the realizations of the random potential, if necessary. In the extreme 
localized case, only one site participates and Sl = 0. In the contrasting limit of fully extended state S'l = 2 — 2/N -^ 2, when 
A^ ^ oo. The participation ratio, hence Sl, exhibits scaling at the metal-insulator transition: P*^^^ = L^^g±{L^/'^w), where 
g± is another universal function. As shown in Fig.|3l the fractal dimension x = 1.4, consistent with the results known in the 
literature. 
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FIG. 2: The minimum of Pb corresponds to the best data collapse. From the figure, we see that y = 2.6 and v — 1.35 are good choices, 
although the minimum with respect to v is very shallow. 
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FIG. 3: The data collapse, as in Fig.Q of the participation ratio of the d = 3 Anderson localization problem. The exponent v and Wc were set 
to the known values 1.35 and 16.5, respectively. The exponent x was then determined to be 1.4. 



The general model of a two-level system (a qubit) coupled to a dissipative heat bath is relevant to a number of interesting 
physical problems. [22] We will consider a specific formulation, the spin-boson model, \TS\ in which the two-level system is 
represented by a spin- 1/2 degree of freedom and the bath is a collection of (bosonic) harmonic oscillators. The Hamiltonian is 



-f^sb — ^-^^^x 



H^ 



^ CTz ^ A„ 



(5) 



where ctx and CTz are Pauli matrices and Hose represents the Hamiltonian of an infinite number of harmonic oscillators. The last 
term couples the z-component of the spin to the coordinates {xn] of the oscillators. The coupling constants {A,i}, together with 
the masses {rrin} and frequencies {tj„}, determine the spectral density of the heat bath, which is given by 



^H = ?E' 



\l 



■ S{U! - LUn) 



mnOJn 

For an Ohmic bath, we can take J{uj) = 2Trauj for uj < ujc and J{uj) = for ui > ujc 



(6) 



At zero temperature, this model has a quantum critical line separating a broken-symmetry phase with (cr^) = Mq ^ from a 
disordered phase with (cr^) = O(seeFig.0J. t^ As with any broken-symmetry phase, its definition requires some care: the states 
with the order parameters +Mq and —Mq are orthogonal and their Hilbert spaces are unitarily inequivalent. An infinitesimal 
external field can pin one of these degenerate vacuua. Thus, the broken-symmetry state has an effective classical description in 
which the environmental degrees of freedom are relaxed around a specific state of the qubit. The corresponding uncertainty is 
zero and so is the von Neumann entropy. In the disordered state, by contrast, infinitesimal field has only infinitesimal effect. 
Right at the quantum critical point, the qubit is maximally entangled with the environment, as will be shown below. So the 
phase transition can be aptly described as a classical-to-quantum transition. There are a number of recent interesting proposals 
to observe this dissipative phase transition using a quantum dot. 124 12511 The renormalization group flows for this transition are 
the same as those of the classical inverse-square Ising model in one dimension; {A/2ujc) plays the role of the fugacity y of the 
kinks in the Ising model, while a maps onto the inverse temperature. Although the linearized renormalization group equations 
are identical to the Kosterlitz-Thouless (KT) transition of the two-dimensional XY model, the physics of the ordered phase is 
entirely different because of the existence of a local order parameter. 
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FIG. 4; The T = phase diagram of a two-level system coupled to an Ohmic heat bath and the corresponding renormalization group flows in 
the regime |(1 — a)| ^ 1 and y <^ 1. The renormalization group trajectories point along increasing values of the flow parameter Z = \n{TLUc), 
where r is imaginary time. The separatrix 2y — —(1 — a) separates a broken symmetry phase ({(Jz) 7^ 0) from the disordered phase 
({o'z} ~ 0). At this critical line the magnetization is discontinuous even though the correlation time (imaginary) diverges. Al a — 1/2 the 
dynamics of the two-level system changes from an overdamped to an underdamped state, but there is no thermodynamic phase transition at 
this point. 



In the ground state of Eq.|5] the reduced density matrix of the spin degree of freedom is determined by the expectation values 

(ax) and (ct^): 



1 









(7) 



After diagonalizing pA, we can easily compute the ground state von Neumann entropy. Hence the behavior of the entanglement 
at the QPT follows directly from the behavior of (cr^: ) and (ctz ) . Since the order parameter (ct^ ) is discontinuous at the transition, 
the von Neumann entropy also jumps by an amount A^. As shown in the Supplementary Information, the magnitude of this 
jump is, to leading order in y. 



AS' = ln2 + (y/2)lny 



(8) 



In the limit of vanishing y, the system goes from being unentangled {{ax) = 0, (az) = 1) to being maximally entangled 
((fa;) = (cz) — 0) as it enters the disordered state. Note that this result depends crucially on a proper treatment of the broken 
symmetry — without the jump in the order parameter, S would be continuous through the transition. We also emphasize that the 
von Neumann entropy is discontinuous (with a singular derivative) even though the correlation length in imaginary time diverges 
at the transition. Because this divergence takes the form of an essential singularity, it does not leave a strong signature in other 
quantities at the critical point. 

In addition to the QPT discussed above, the spin-boson model also undergoes a dynamic crossover at a = 1/2, below which 
the spin can exhibit damped coherent oscillations on short time scales. For a > 1/2 the dynamics of the disordered phase is 
characterized by an overdamped decay. If the central claim of this paper is correct, the von Neumann entropy should be analytic 



at a = 1/2 because no phase transition occurs at this point. This appears to be consistent with numerical renormalization group 
calculations of the anisotropic Kondo problem. 12611 

We now turn to a simpler but equally important example, the damped harmonic oscillator Consider a single harmonic 
oscillator (momentum p, position x, mass M, and frequency ujq) coupled to an environment that also consists of harmonic 
oscillators. Again we examine the case of Ohmic dissipation, where the spectral function (defined as in Eq.|6j is J{u!) = rjuj for 
small uj and zero above the cutoff tJc- The ground state expectation values of x^ and p^ are given by 12711 

where k — ri/2Mu!o is the friction coefficient and 



^^ (9) 

1^0 ' 



m^ / lnr + f_l (10) 

TTV K^ — 1 V K — VK^ — 1 / 



At K = 1 the system crosses over from damped oscillatiory to overdamped behavior (Vk^ — 1 is to be replaced by iVl — k^), 
just as the spin-boson model does at a = 1/2. The function /(k) is real for all k > and has identical power series expansion 
whether we approach k = 1 from the positive or the negative side. Thus it is analytic at k = 1. 

At zero temperature, the normalized reduced density matrix for the damped harmonic oscillator has matrix elements 12711 

{x'\pa\x") - v/46Ae-"("'-"")'-''("'+"")' (11) 

where a — {p^)/2h^ and b — l/8(a;^). To compute the von Neumann entropy, we first note 12811 that 

Tr{pA\npA)= lim ^ f dx' {x'\p\\x') (12) 

n^i on J 

The subsequent derivation of S is reproduced in the Supplementary Information; the result is 

2 \a-b) 2\ b \^-VbJ 

Note that S* — > as 6 -^ a, which corresponds to the minimum uncertainty \J {x'^){p'^) ~ h/2. The uncertainty relation is 
satisfied only for b < a. We have verified that S is analytic at k = 1 (see Supplementary Information) — this follows from the 
analyticity of the function /(k). 

Despite the popularity of field-theoretic perturbative methods involving Feynman diagrams, some of the major breakthroughs 
in condensed matter physics have come from formulating wave functions that reflect the unique correlations of a quantum 
system with infinitely many degrees of freedom. The BCS wave function is remarkable because it incorporates the broken 
global gauge symmetry. The Bijl-Dingle-Jastrow-Feynman wave function for liquid ^He contains the essential ingredients of 
interacting superfluid matter. Similarly, the Laughlin wave function explains not only the fractional quantum Hall effect, but 
also the fractionalization of charge. We believe that the non-analyticity of the von Neumann entropy is a useful criterion for all 
quantum critical points at which the character of the ground state wave function changes. It provides a unique insight into the 
internal correlations of a many -body wave function. The approach described in the present paper can be applied to a wide variety 
of QPTs for which the conventional energy criterion fails. In particular, it can be immediately applied to transitions between 
integer quantum Hall plateaus. ll29ll 

SUPPLEMENTARY EQUATIONS, FIGURES AND NOTES 
Derivation of AS for the spin-boson model 

To compute AS*, we need to know the expectation values (ax) and (ct^) at criticality. Adapting the results of Anderson and 
Yuval 13011 . we find that (ctz) jumps by a/T/c^ along the quantum critical line 2y = — (1 — a). The expectation value of cTx is 
most easily obtained by differentiating the ground state energy: (ax) = — 2 (dE/dA). Since the ground state energy maps onto 
the free energy of the inverse-square Ising model, we can make use of the free energy scaling relation Ii30n31il 

dF{l)^y^{l)e~^dl (14) 



Here F is the free energy in units of Wc, so we have the direct correspondence (a^) — —2 (dE/dA) -^ —dF/dy. 

On the separatrix 2y = —(1 — a), the flow of the fugacity is given by y{l) = y/{l + 2yl). Since y — > as Z ^ oo, 
the renormahzation group equations (vaHd for y ^ 1) become asymptotically exact in this limit. The fixed point at ^ = cx3 
corresponds to the fully polarized state of the Ising chain, for which the free energy is zero. So if we tune to criticality, we can 
obtain a closed-form solution for F by integrating Eq.[2]all the way to / = cx). Differentiating with respect to y then gives 



(^.> = I 



1- — 



1 



=i/2y 



2y {^yy 



dt- 



1/2!/ 



(15) 



In the limit y ^ 0, Eq.llSlassumes the form (ctj.) w 2y, while the discontinuity in (cr^) goes as 1 — y. Using the reduced density 
matrix given in the text, it is easy to show that 



A5 = ln2+(y/2) hi y + 0(y) 



(16) 



Derivation of S for tlie damped harmonic oscillator 

As described in the text, the von Neumann entropy of the damped harmonic oscillator is given by 



S = -TY{pA\npA) = - lim — / dx'{x'\p\\x') 
ri^i on J 



(17) 



After inserting n—\ complete sets of position eigenstates, we can use the expression for the matrix elements to recast the integral 

as 



dx'{x'\pWx') = 



— j dxi---dxn e-"'*^--^ 



(18) 



M is a tri-diagonal matrix with Ma — 2{a + b) and A/^+i^i = Mi^i+i — —{a ~ b), and therefore has eigenvalues 7„ 
2(a + b) — 2(a — b) cos(27rm/n), with ni — 1, . . . ,n. So the von Neumann entropy becomes 



lim 



d 



*i dn 



{^br" n > 



-1/2- 



For convenience, we write 



J]7,„ = 2"(a-fe)"P 



(19) 



(20) 



where the product P is defined as 



^-n 



27r7Tl\ 

a — cos I 1 



a + 6 
a — b 



Our main task is to evaluate this product by performing a suitable integral transform. Note that 



dlnP 
da 



E 



1 



1 



dz ■ 



1 



nz 



n-l 



;a-cos(2ip) 27nTc a-\{z + \) z^-l 



(21) 



(22) 



provided the contour C encloses the poles at Zm 

a In P 2n 



^i2-Kinjn 



z± — a± \/ o? — 1, we obtain 



(Fig-E). Deforming the contour to pick up the additional poles at 



da 



z+ — z- \ zV: — 1 z" — 1 



2n z!f + 1 



z+ — z_ z? — 1 



(23) 



where the second equality follows from z+2_ — 1. Integrating this result over a yields a closed-form expression for P. The 
integral is most easily evaluated by transforming to the variable z_|_: 



InP 



dzu 



n zl + 1 



Z I z" - 1 



In ( z« + — - 2 ) + const. 



(24) 




FIG. 5: The contour C encloses the poles on the unit circle at z^ = g'^jrm/n -^jjjj a > \ (h < a), the additional poles at 2;± = a ± \J d^ — 1 
are located on the real axis. 
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FIG. 6: Von Neumann entropy (S) of the damped harmonic oscillator with u)c ~ 100 a;o. Both S and its first derivative (inset) are perfectly 
smooth through the dynamic crossover at k = 1. Higher-order derivatives are similarly well-behaved due to the analyticity of /(k) at k = 1. 



This is consistent with Eq.|^if we set the integration constant equal to — n In 2. In terms of the original parameters a and b, the 
product is given by 



1 
P= — 

2" 



Vb' 



Vb' 



/a—VbJ \\/a + VbJ 

and the rest of the calculation follows trivially from Eqs.^]and|20| The von Neumann entropy is 



(25) 



2 \a-b 



1 la 



In 



Vb' 



\/a- \/6 



(26) 



where the ratio b/a depends on the friction coefficient n and the ultraviolet cutoff cjc (recall that the uncertainty principle requires 
b < a). S is an analytic function of k as long as b/a is analytic — which, in turn, is guaranteed by the analyticity of the function 
/(k). The plots in Fig.|6lprovide further proof of this point. 
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